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TECHNICAL MEMORANDUM 


THEORETICAL CONSOLIDATION OF ACOUSTIC DISSIPATION 

1. INTRODUCTION 


In many engineering problems, the effects of dissipation can be extremely important. Dis- 
sipation can be represented by several parameters depending on the context and the models that 
are used. Some examples of dissipation-related parameters are damping ratio, viscosity, resis- 
tance, absorption coefficients, pressure drop, or damping rate. This Technical Memorandum (TM) 
describes the theoretical consolidation of the classic absorption coefficients with several other dis- 
sipation parameters including linearized resistance. The primary goal of this TM is to theoretically 
consolidate the linearized resistance with the absorption coefficient. As a secondary goal, other 
dissipation relationships are presented. Table 2 in section 6 presents all the dissipation relationships 
in this TM. 


Reducing the linearized resistance parameter as a function of acoustic absorption coeffi- 
cient is a primary goal because of its importance in impedance or lumped parameter modeling of 
the dynamics in pipe systems without flow. Examples of its use would be for dead-headed cavities, 
closed pipes, or sensor ports when the flow resistance is reduced to zero. It is well known that pres- 
sure drop is a function of flow rate, and linearized resistance follows directly from a Taylor series 
expansion of Ap(m) as shown in equation (1): 


A p(m) ■ 


Ap(rii) + r)A,> 


dm 


[m - m^ + Oyn' 2 ^ 


( 1 ) 


The dependent variables in this TM are made up of a steady part and an oscillatory part; e.g., 
p = p + p, so in oscillatory form, equation (1) can be written as equation (2) after dropping 

the higher order terms, 0\m 2 ) : 


Ap(m)' = 


dA p 

dm 


m 


(2) 


Linearized resistance can be defined as the linear relationship between oscillatory pressure 
drop and mass flow rate and is shown in equation (3): 




dAp 

dm 


( 3 ) 


In a system without flow, the usually dominant form of resistance due to steady flow viscous 
effects, orifice pressure drops, or other dissipative components, is equal to zero. To model the effects 
of dissipation in a fluid system without flow, acoustic absorption coefficients can be used. The 
absorption coefficients, described in further detail in section 4, are used classically as estimates of 
wave propagation dissipation in a fluid. They are useful because they can be estimated directly with 
known properties of the fluid and also can be related to a linearized resistance so they can be used 
in a dynamic system without flow. 

A resistor, inductor, and capacitor (RLC) circuit is commonly used as an analogy in fluid 
systems. There are several RLC fluid analogies, however one of the most commonly used analogy 
considers the RL in series with a shunt C for each component or lump. Reference 1 (chp. 12), for 
example, in dynamic fluid systems analysis, uses this formulation. Reference 2 (pp. 270-275) also 
discusses the circuit analogy in terms of modeling the feed system components in the combustion 
stability problems. Reference 3 (pp. 18-22) describes this fluid analogy along with others including 
the forms of their equations. Recognizing the dissipation formulation and wave equation used is an 
important step in consolidating dissipation parameters. 
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2. FORMULATION DESCRIPTION 


The key to consolidating the dissipation terms is to develop identical theoretical formula- 
tions while incorporating generic forms of the dissipation parameters. By casting each theoretical 
formulation as a damped wave equation, each dissipation parameter can be directly compared. 
Keeping the dissipation terms generic from the start of the theoretical development, rather than 
substituting in known parameters, provides a good approach. In this section, a damped wave equa- 
tion is developed with linearized resistance and one is also provided with the acoustic absorption 
coefficient. 


A force balance in a differential section of pipe can be simplified to equation (4). This is 
similar to that described in reference 1 (p. 21, equation (2.1)) except that the pipe slope is ignored. 
The significant observation in this formulation is that the dissipation is represented by shear stress 
rather than viscous terms as in the Navier-Stokes equation: 


A d -P + Bd. 


du dll 

— h U 


dz g c y dt dz 


+ TP = 0 


(4) 


As described in reference 1 (p. 21), the usage of shear stress, a steady flow resistive force, is an 
assumption typically used in transient flow calculations and forms the basis for the unsteady flow 
equations. Also described in reference 1 (p. 290) and reference 4 (p. 231), it is recognized that the 
wall shear stress is not in phase with the mean velocity in pulsatile flow and the inertia term must 
be modified. Additionally, using wall friction as a static function of the mean velocity underesti- 
mates the wave attenuation at moderate and high frequencies (ref. 1, p. 290). Considering that the 
advective acceleration term, uxdu/dz, is small in low Mach number unsteady flows, it is eliminated 
from the equation; equation (4) becomes the following momentum equation: 


. dp pA du n A 

A- + - + Pz = 0 

dz g c dt 


(5) 


In steady flow, a force balance gives equation (6), except now it is assumed that this relationship 
also applies in transient flow, thus the dependent variables contain an oscillatory component: 

A pA = tLP . 


In steady state, equation (6) reduces to A pA = z PL. Therefore, equation (6) can be reduced to 
oscillatory form as equation (7): 


A p A = z'L P . 


(7) 
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In oscillatory form, equation (5) can be written as equation (8): 


A ¥ + pAM pz , = () 

dz g c dt 


( 8 ) 


The relationship, A p = R L m, holds from equation (2). After substituting the oscillatory mass 
flow rate, m = pAu, into this relationship, the oscillatory pressure drop, A p, can be substituted 
into equation (7). The result is shown as equation (9) and is the oscillatory shear stress. Recall the 
hydraulic diameter, D H = 4 A/P. 


, _ pAD H R L 

W 


(9) 


Equation (9) can be substituted directly into equation (8). This is the final form of the oscillatory 
momentum equation while keeping the linearized resistance as a generic parameter: 


dp' p du pARj , „ 

— + — + - —u = 0 . 

dz g c dt L 


(10) 


The approach used by reference 1 is to obtain shear stress specifically from the Darcy- 
Weisbach equation: 


&p = f 


T „ 2 

L pu 
d h 2 Sc 


( 11 ) 


And then by substituting equation (11) into equation (6), the shear stress relationship in equa- 
tion (12) is obtained. Since the shear stress will be used in oscillatory flow, the absolute value 
on the velocity term ensures that the shear stress will always oppose the direction of velocity. 


T = 


pMu 
8 g c 


( 12 ) 


Equation (12) is then linearized to obtain equation (13): 


T' = f^-U' 

4g c 


(13) 


At this point, equation (13) is used directly in equation (8) to obtain a final form of the oscillatory 
momentum equation, which does not contain the generic linearized resistance parameter as shown 
in equation (14). This form of the oscillatory momentum equation is not helpful for determining 
linearized resistance since it has already been eliminated in the equation. 
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As an aside, equation (11) with equation (3) can also be used to obtain the linear resistance in 
a pipe with flow shown as equation (15): 


R L = f 


L m 
D n pA 2 g c 


( 15 ) 


The next step in obtaining a damped wave equation with a general linearized resistance term 
is to bring in the continuity equation. Also considering that the advective term is small, the oscil- 
latory continuity equation can be written simply as equation (16) and is also limited to low Mach 
number flows: 


dz pc/Q dt 


(16) 


By differentiating equation (10) with respect to z, differentiating equation (16) with respect to t, 
and combining to eliminate the velocity terms, the following wave equation is obtained. This is 
a damped wave equation that incorporates the linearized resistance: 

<jV 1 d 2 p AR L g c dp' = 

dz 2 ciq dt 2 %L dt 


Reference 5 provides the development of a damped wave equation for waves in a viscous 
stationary medium. The linear momentum equation is derived by incorporating the pressure drop 
per unit length due to visco-thermal friction and can be written as equation (18): 


dp' p du 2 paan . n 

— + — +— — -u = 0 

dz g c dt g c 


(18) 


And similarly combining with continuity, the damped wave equation becomes equation (19): 


d 2 p' 1 d 2 p' 2 a dp' 

dz 2 dt 2 % dt 


(19) 


By direct comparison of the damped wave equations, equations (17) and (19), the relationship 
between linear resistance and absorption coefficient can be deduced and is described in section 3. 


3. DISSIPATION PARAMETER CONSOLIDATION 


Several dissipation parameters are consolidated in this section. The linearized resistance 
(using pressure and mass flow rate), pressure damping rate, energy damping rate, damping ratio, 
and linearized resistance (using head and flow rate) are all related to the acoustic absorption coef- 
ficients. The relationship between each individual parameter can also be found by comparison 
of the parameters in table 2 in section 6. 

3.1 Linearized Resistance (Pressure and Mass Flow Rate) 

By comparing equations (17) and (19), an equality can be written as equation (20): 

ar l 8 c = 2a ( 2 i 

OqL a 0 

Solving for the linear resistance gives the following relationship: 



2 LaciQ 

A S C 


( 21 ) 


Equation (21) gives the linearized resistance as a function of the absorption coefficient. Absorption 
coefficient is in the form of dissipation per unit length with units of neper per unit distance. The 
absorption coefficient can be evaluated as described in detail in the next section. 


3.2 Pressure Damping Rate 

Another common parameter used for dissipation is the damping rate. The damping rate can 
be in two forms and care must be taken when using this parameter. The energy damping rate 
is twice that of the pressure damping rate because of the relationship between total energy and 
pressure. This wifi be described further below. 

Because this derivation is usually discussed only in part in literature, 4 ’ 6 ’ 7 it is shown here 
in completion. To represent damping, the frequency or the wavenumber can be regarded as com- 
plex. For damping rate, the frequency is regarded as complex where co — xb-kt^. In this form, 
o p > 0 indicates that the wave decays as is seen in equation (22): 

pit) = ) = p Qe -tot e i (22; 

The acoustic energy density of a plane wave can be represented as equation (23): 
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(23) 


e(t) = 


P(t) 1 g c 

pal 


Since the complex form of equation (22) is normally used as a convenience in representing 
a harmonic signal in acoustics, the complex conjugate must be incorporated in the evaluation 
of p(t) 2 . Also assuming that the fractional change of amplitude is small in one cycle of oscilla- 
tion, | o p / c»\ « 1 , the quantity e ° pt can be regarded as constant in time. The evaluation of pit) 2 
becomes equation (24): 


pitY 


-Re 

2 


r 


D 2 ' 
P 0 e 


2i(a>t-io.,t-kz ) 2 -2a„t 

v y ’+Pde 1 


V 


(24) 


And the energy density can be simplified to equation (25): 


e(t) = 


d 2 -2 at 

P 0 e p g c 


2pa {) 


(1 + cos(2 cot - 2 kz)) . 


(25) 


Continuing with the small fractional change in amplitude assumption, the quantity e ° pt can 
again be regarded as constant in time and factored out of the integral in equation (26). The period 
is over one cycle and can be written as T- 2nla>. Parseval’s theorem can also be applied directly. 


£ = — 
T 


1 1 

pmdtPj 


pjtTgc 
P a 0 j 


dt = 


Pte 


-lot 


Sr 


2p«o 


It can be shown that the following equality is satisfied by evaluating: 


1 d(e) 
2 (e) dt 


(26) 


(27) 


For the absorption coefficient, the relationship in equation (28) holds. In this case, the wave 
number is regarded as complex where k^k-ia. In this form, a> 0 indicates that the waves decay 
as is seen in equation (28): 


p{t) = P 0 e- io)l e ikz e- az = P {) e~ ia)t e ikz e^ 1 . 
Using equation (28) and following the same procedure gives equation (29): 

1 die) 

a= — — ^ . 

2 a 0 (e) dt 


(28) 


( 29 ) 
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And so the following relationship is obtained by comparing equations (27) and (29): 

o p = aa {) . (30) 

3.3 Energy Damping Rate 

Acoustic energy can be described as the energy density multiplied by volume. Assuming the 
oscillatory part of acoustic energy follows the form it is seen from equation (25) that 

the energy damping rate is given by equation (31): 

o E = 2o p . (31) 

Or in terms of absorption coefficient, it is given by equation (32): 

o E = 2aa 0 . (32) 

3.4 Damping Ratio 

The damping ratio is used at times and is commonly used in finite element modeling because 
of analytic features designed to couple acoustics and structures. It is also representative of the 
damping parameter used in a system with a single degree of freedom. Analogous to a mass-spring- 
damper system, the wave equation can be rearranged and written as equation (33). The parameters 
m, b, and k do not represent the actual mass, spring constant, and damping coefficient. For the 
parameters to be representative, the wave equation would be written as a force balance, however 
doing this would only introduce a linear factor that would be eliminated when rewriting in a form 
for an analogous equation comparison: 


m^- + b^ + kp' = Q. (33) 

dt at 

The definition of damping ratio of a damped harmonic oscillator can be represented by the damp- 
ing coefficient, equation (34): 


b = 2 ^kmt, = 2 mco £ . (34) 

To provide a direct comparison, the lossless wave equation parameters for m and k, in the form 
of equation (33), are taken from equation (19). Also, by incorporating equation (34) using the 
newly defined parameter for m, the damped wave equation with damping ratio can be written as 
equation (35): 


f 


\ 


2 

a 0 J 


~\2 / 

a p 


dt z 


+ 


r 




2co 

«0 


dp_ 

dt 


+ 


dz 7 


0. 


(35) 
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By direct comparison of equations (19) and (35), the equality shown as equation (36) is obtained: 

-2a = -2co £ (36) 

a () CIq 

Damping ratio as a function of absorption coefficient can be related as equation (37): 


£ = . (37) 

(O 

3.5 Linearized Resistance (Head and Flow Rate) 


Care must also be taken for resistance terms in various formulations. For example, it is com- 
mon to normalize resistance to density or length. Additionally, the linearized resistance may be 
based on head and flow relationships. For example, in reference 1 , the linear resistance is based on 
head and flow rate, rather than pressure and mass flow rate and is also normalized to length. Their 
momentum equation is given by equation (38): 


dti 
— + 
dz 


J_d£ 

gA dt 


+ R 


'l,w c J 


0. 


(38) 


With ti = —p’ 
Pg 


and q = Au, equation (38) can be written as equation (39): 


V + £_M + Pg AR LW u , = Q (3 , 

g c gc 

And by comparing the momentum equations, equations (10) and (39), the equality given by equa- 
tion (40) is obtained: 


pAR[ _ pgA R/ jy 

L gc 


(40) 


Thus, the relationship between the Wyle & Streeter 1 form of linearized resistance (using head 
and flow rate) and the linearized resistance (using pressure and mass flow rate) is described by 
equation (41): 


r l,w 


RLgc_ 

L g 


(41) 


It is also clear that the Wyle & Streeter 1 form of linearized resistance is normalized to length. 
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In terms of absorption coefficient, the Wyle & Streeter 1 linearized resistance relationship is given 
as equation (42): 


r l,w 


2aciQ 

gA 


(42) 
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4. ABSORPTION COEFFICIENTS 


Sources of dissipation can be broken into two main categories — those intrinsic to the 
medium and those associated with the boundary of the medium. The losses intrinsic to the medium 
can be further divided into three basic types: (1) viscous losses, (2) heat conduction losses, and 
(3) losses associated with internal molecular processes. At the boundary there are viscous losses 
and losses to heat conduction between the fluid and the surface. Without derivation the absorp- 
tion coefficients are provided in this section. Most acoustics textbooks discuss these absorption 
coefficients in further detail. 3 ’ 4 For the dissipation associated with the boundary, the absorption 
coefficients are usually derived for standing waves. However, reference 7 gives the formulation of 
damping rate that includes nonresonant conditions. The forms in reference 7 are provided in this 
TM to account for both resonant and nonresonant conditions. 

The total absorption coefficient can be regarded as the sum of the absorption coefficients 
for the individual loss mechanisms. Superposition is typically justified in practice, and true when 
losses are small; however, in general, they do have interactions with each other. The unit for the 
absorption coefficient is neper per unit length where a neper is analogous to a bel, but represents 
the natural logarithmic scale rather than the decadic logarithmic scale. 


The losses intrinsic to the medium can be further divided into three basic types: (1) viscous 
losses, (2) heat conduction losses, and (3) losses associated with internal molecular processes. The 
complete absorption coefficient is defined as equation (43), but the combination of equations (44) 
and (45) under Stokes’ assumption is regarded as the classical absorption coefficient. 

4.1.1 Intrinsic Absorption Due to Viscosity 

The intrinsic absorption in a fluid due to viscosity is given by the intrinsic viscous absorp- 
tion coefficient, equation (44). In Stokes’ assumption, the bulk volume viscosity is zero. Substantial 
theoretical development in the field has occurred since Stokes’ pioneering work that has introduced 
the bulk volume viscosity parameter into the absorption formulation. The sound speed dispersion 
for intrinsic viscous effects is very small and usually too small to measure. 


a = ^ a, = + a K + a M + a 



(43) 


4.1 Intrinsic Dissipation 



(44) 
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4.1.2 Intrinsic Absorption Due to Thermal Conduction 


The intrinsic absorption in a fluid due to thermal conduction is given by the intrinsic ther- 
mal absorption coefficient, equation (45). Variations in temperature occur between the slightly 
warmer compressed portion and slightly colder expanded portion of the wave. The variation in 
temperatures gives rise to heat flow as the fluid medium tries to reestablish temperature equilib- 
rium. The heat flow reduces the amount of energy available for the sound wave and represents an 
acoustic loss. The wave due to this absorption mechanism is shown to be nondispersive. For most 
gases the absorption due to intrinsic viscosity and intrinsic thermal conduction is comparable; for 
most liquids the intrinsic absorption due to thermal conduction is usually negligible. 


co 2 (y~1)k 


2pciQ Cp 


(45) 


4.1.3 Intrinsic Absorption Due to Molecular Thermal Relaxation 

The intrinsic absorption in a fluid due to molecular thermal relaxation takes into account 
the internal structure of the molecules and the interactions between them that lead to internal 
vibrations, rotations, ionizations, and short-range ordering. In this case, the pressure is not just 
dependent on the local instantaneous values of density and pressure; additionally, it has a depen- 
dence on the rate at which density and temperature are changing. The molecular thermal relaxation 
time is associated with the time needed for a new equilibrium to be established. Transformation 
of the mechanical energy into heat represents a loss to the sound wave. The absorption is most 
significant at the molecular relaxation frequencies of the fluid. The fluid may have several relax- 
ation frequencies in which the total molecular thermal relaxation is the sum of all the absorption 
coefficients calculated separately for each of the relaxations. The sound speed transitions from 
the equilibrium sound speed, a () , at low frequencies that are well below the relaxation frequency, 
to the frozen sound speed, at high frequencies that are well above the relaxation frequency. 
Experimental data for a particular fluid are needed to obtain the parameters, /r max and f M , which 
are the maximum value of absorption per wavelength and the relaxation frequency, respectively. 
These parameters are both functions of temperature. This form of absorption can be responsible 
for most of the intrinsic sound absorption in the audio and low frequency regions. The absorption 
coefficient due to molecular relaxation is given by equation (46) and the phase speed dispersion for 
molecular relaxation is given as equation (47): 


and 


a M ~ 


2fi 


max 


«0 


1 + 


Jm 

Jm 

f 


\2 


a 


’ph a 0 


1 + 


( 2 r\ 2 

a oof 

%Jm 


1 + 


f r \2- 

«oo / 


(46) 


(47) 
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v a 0 Jm J 


4.2 Boundary Dissipation 


Boundary dissipation is an entirely separate mode of absorption due to the passage of 
a wave over a surface boundary. The surface exerts a frictional shear force on the overlying fluid 
while thermal conduction allows heat transfer to take place between the fluid and the surface. Both 
result in loss of energy from the wave. 

Similar to a boundary layer thickness in steady flow, acoustic waves have an analogous 
acoustic boundary layer thickness. It is much thinner because the boundary layer only has a period 
of half a cycle to develop. For the formulas to hold, the acoustic wavelength must be much greater 
than the acoustic boundary layer thickness A » <5 . The viscous and thermal acoustic boundary 
layer thicknesses are given as equations (48) and (49), respectively: 



and 


<S, 



(48) 


(49) 


For the wavelength, A = 2na () lo), to be much greater than the acoustic boundary layer thickness, 
equations (50) and (51) must be satisfied: 


and 


i 

7 TCIq Y 2 p 


(50) 


1 I j-10) 
7ia {) ]]2pPr 


(51) 


The combination of the wall viscous absorption coefficient, equation (52), and the wall ther- 
mal conduction absorption coefficient, equation (54), is regarded as the combined absorption 
coefficient for wall losses. Normally, the formulation of the following boundary layer absorption 
coefficients are derived for resonant conditions, however, reference 7 gives the formulation of damp- 
ing rate that includes nonresonant conditions. The absorption coefficients at general conditions are 
provided so that both resonant and nonresonant conditions can be estimated. For other pipe shapes, 
the radius can be replaced with a hydraulic diameter term, R=D H I2. Reference 8 compares experi- 
mental data, to a very fine degree, to the combined absorption coefficient theoretical formulation. 

4.2.1 Viscous Losses at a Rigid Wall 

The boundary layer absorption in a fluid due to viscosity for a cylindrical pipe is given by 
the wall viscous absorption coefficient, equation (52). This absorption coefficient accounts also 
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for nonresonant conditions. The phase speed dispersion for wall viscous effects is given as equa- 
tion (53). It shows that the sound speed is slightly slower. 


a , 


1 Ijuco f sin(2/cL) ' 


( 52 ) 


a 0 R]j 2p { 2kL 


and 


a P h ~ a 0 1 

y co ) 


1 


(53) 


4.2.2 Thermal Conduction at an Isothermal Wall 

The boundary layer absorption in a fluid due to thermal conduction for a cylindrical pipe is 
given by the wall thermal absorption coefficient, equation (54). This absorption coefficient accounts 
also for nonresonant conditions, but also end-wall effects that are important in thermal conduc- 
tion. The sound speed dispersion for wall thermal effects is given as equation (55). It shows that 
the sound speed is slightly slower. 


1 


pcof y-1 V sin(2 kL) R N 


— , 1 H 1- It — 

V VPr ) \ 2 kL L ) 


(54) 


a 0 R ]/ 2 p 


and 



(55) 
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5. EXAMPLE 


As an example, a 5-in-long senseline with an inner diameter of 0.5 in is examined. A sche- 
matic is shown in figure 1 . The contents are air at a temperature of 540 °R and pressure of 4 torr. 


D= 0.5 in 


I 



L = 5 in 


Figure 1 . Sensor port example. 


An expression for complex pressure ratio, p^/p^, is needed where A is the open end and B 
is the closed end. From reference 1, for a simple pipe closed at one end, this complex pressure ratio 
is described with the hyperbolic secant as equation (56): 

— = sech (yL) . (56) 

Pa 


And the propagation constant, 7, can be written as equation (57): 


7 = 


r 


co 

a 0 ) 


^ ( (DgAR LW ^ 


+ 1 


2 

a 0 


(57) 


In this case, the acoustic dissipation is dominated by the viscous and thermal losses at the senseline 
wall. Since the wall effects dominate in this problem, the acoustic absorption coefficient is found by 
substituting the fluid parameters in table 1 and the geometric parameters into equations (43), (52), 
and (54). 
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Table 1. Fluid parameters for air with a temperature of 540 °R (80.33 °F) 
and pressure of 4 torr (0.077347 psia). 


Parameters 

Variable 

Value 

Units 

Viscosity (dynamic) 

P 

1.0373x10-6 

lb f /in-s 

Density 

P 

2.2367x10 -7 

lb m /in 3 

Sound speed 

a 0 

13,671.6 

in/s 

Specific heat ratio 

r 

1.4 

- 

Prandtl number 

Pr 

0.70641 

- 


Additionally, substituting these parameters into equation (42) to obtain the linearized resis- 
tance term, using equation (57) to obtain the propagation constant, and then lastly combining into 
equation (56), the complex pressure transfer function can be obtained. Equation (42) allows access 
of the absorption coefficient in the form of linearized resistance and thus can be applied directly 
into the complex pressure ratio model. Figure 2 shows both acoustic absorption coefficient and 
damping ratio as a function of frequency. Figure 3 shows the gain and phase of the complex pres- 
sure transfer function between location A and B due to acoustic dissipation. 



0.25 


- 0.2 


O 

O) 

3 


73 

Q) 


o 


-0.05 


Frequency (Hz) 

Figure 2. Acoustic dissipation parameters. 
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Frequency (Hz) 



Frequency (Hz) 


Figure 3. Pressure transfer function for a sensor port using acoustic absorption: 
(a) Gain of p B lp A and (b) phase of p B ! P a- 
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6. SUMMARY 


This TM describes the theoretical consolidation of the classic absorption coefficients with 
several dissipation parameters including linearized resistance. Table 2 provides the various relation- 
ships developed in this TM in terms of absorption coefficient. This allows engineering problems 
that use various dissipation parameters to account for acoustic dissipation by easily estimating and 
incorporating the needed parameters into their particular models. This TM summarizes the con- 
solidation of several dissipation parameters and gives a simple example of a senseline problem. 


Table 2. Acoustic dissipation relationship as a function of absorption coefficient. 


Linearized resistance (pressure/mass flow rate) 
R l =^ 

m=m 

2 Laa^ 

R i — 

a 8c 

Pressure damping rate 

o p = aciQ 

Energy damping rate 

o £ — Iccciq 

Damping ratio 

ccciq 

CO 

Linearized resistance per unit length (head/flow rate) 

1 dAh 

r l,w = ~ 

L da 

1 9 = 9 

2 aa 0 
k l,w - , 

g A 
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